DRAFT VERSION FEBRUARY 5, 2008, ACCEPTED FOR PUBLICATION IN THE ASTROPHYSICAL JOURNAL 
Preprint typeset using L^T^X style cmulatcapj v. 12/14/05 



A SEARCH FOR H 2 IN THE STRONGLY LENSED QSO MG 0751+2716 AT Z = 3.2 

DOMINIK A. RlECHERS 1 , AXEL WEISS 2 , FABIAN WALTER 1 , CHRISTOPHER L. CARILLI 3 , AND KlRSTEN K. KNUDSEN 1 
draft version February 5, 2008, accepted for publication in the Astrophysical Journal 

ABSTRACT 

We present a search for 183 GHz H 2 0(3i3^2 2 o) emission in the infrared-luminous quasar 
MG 0751+2716 with the NRAO Very Large Array (VLA). At z = 3.200 ± 0.001, this water emis- 
sion feature is redshifted to 43.6 GHz. As opposed to the faint rotational transitions of HCN (the 
standard high-density tracer at high— z), H20(3i3-^22o) is observed with high maser amplification 
factors in Galactic star-forming regions. It therefore holds the potential to trace high-density star- 
forming regions in the distant universe. If indeed all star-forming regions in massively star-forming 
galaxies at z > 3 have similar physical properties as e.g. the Orion or W49N molecular cloud cores, 
the flux ratio between the maser-amplified H20(3i3^22o) and the thermally excited 12 CO( J=1^0) 
transitions may be as high as factor of 20 (but has to be corrected by their relative filling fac- 
tor). MG 0751+2716 is a strong 12 CO(J=4^3) emitter, and therefore one of the most suitable 
targets to search for H20(3i3— >22o) at cosmological redshifts. Our search resulted in an upper 
limit in line luminosity of ijj 2 < 0.6 x 10 9 K kms _1 pc 2 . Assuming a brightness temperature of 
Tb(H20) ~ 500 K for the maser emission and CO properties from the literature, this translates to a 
H20(3i3— >2 2 o)/ 12 CO(J=4^3) area filling factor of less than 1%. However, this limit is not valid if 
the H20(3i3-^22o) maser emission is quenched, i.e. if the line is only thermally excited. We conclude 
that, if our results were to hold for other high-z sources, H2O does not appear to be a more luminous 
alternative to HCN to detect high-density gas in star-forming environments at high redshift. 
Subject headings: galaxies: active, starburst, formation, high redshift — cosmology: observations — 
masers 



1. INTRODUCTION 

Over the past years, massive amounts of dust and gas 
have been detected in distant quasars, allowing to study 
the properties of molecular gas in the early epoch of 
galaxy formation. Molecular gas in high-redshift galax- 
ies is commonly traced by CO emission and was found 
in >30 galaxies at z > 2 to date, out to the highest red- 
shift quasar , SDSS J1148+ 5251 at z = 6.42 (Walter et al. 
120031 120041 Bertoldi et al. l2003|) . The observed molecu- 
lar gas masses in excess of 10 10 M Q provide the requisite 
mate rial fo r star formation (SF, e.g. Solomon & Vanden 
Bout l2f)05l . 

The presence of abundant molecular gas, the fuel for 
SF, has led to the hypothesis that the tremendous far- 
infrared (FIR) luminosities (> 10 12 L Q ) of these high- 
redshift objects are not only powered by active galactic 
nuclei (AGN) but also by major starbursts (SBs), sug- 
gesting that this population represents the formation of 
large spheroidal galaxies beyond redshift 2 (Blain et al. 
2002). However, the relative contribution to the FIR lu- 
minosity from dust heated by SF and AGN activity in 
high-z sources remains subject to discussion (e.g. An- 
dreani et al. l2003|) . 

While the lower- J transitions of CO are a good indica- 
tor for the total molecu lar gas conte nt of a system (e.g. 
Carilli et al. 119991 12003 Riechers et al. l2006|) . they can be 
excited at relatively low densities; the critical density lies 
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at only nn 2 — 10 3 cm~ 3 . Hence, it is a relatively poor 
tracer of the denser gas directly involved in massive star 
formation. The standard tracer of the dense molecular 
gas phase is HCN, and the critical density to excite its 
lower-J transitions is 71h 2 — 10 5 cm -3 (due to a higher 
dipole moment in comparison to CO). In this context, 
recent studies of the dense molecular gas phase in local 
(z < 0.3) luminous and ultra-luminous infrared galaxies 
(LIRGs/ULIRGs) have shown a correlation between the 
HCN luminosity and the star-formation rate (SFR) as 
traced by the FIR luminosity. This correlation is much 
tighter than th at betw een the CO and FIR lum inositie s 
(Solomon et al. l1992at Gao & Solomon 12004^1 I2004h1. 
Thus, observations of HCN towards high- z sources hold 
the potential to pin down the contribution of SF to the 
total FIR luminosity. Unfortunately, emission lines from 
HCN in LIRGs/ULIRGs are typically by a factor of 4 
to 10 fainter than those of CO (even a facto r of 25 - 
40 in ordinary spiral galaxies, Gao & Solomon [2004b) , 
which hinders systematic HCN surveys at cosmological 
distances using current telescopes. Until today, deep ob- 
servations of HCN at z > 2 resulted in four detections 
and four upper limits (Solomon et al. 2003: Vanden Bout 
et al. I2004T Carilli et al. l2005t Wagg etal. 120051) . 

Due to the faintness of emission connected with ther- 
mally excited rotational transitions of HCN, it is im- 
portant to investigate whether another physical process 
can be found which produces emission lines with signifi- 
cantly higher luminosities. The 183 GHz 3i3 — > 220 emis- 
sion line of para-H20 holds the potential to be such a 
tracer: in warm, dense SF regions (riH 2 > 10 5 cm~ 3 ), 
this line is collisionally pumped at relatively low ki- 
netic temperatures (Tkin = 50 — 100 K, e.g. Cernicharo 
In contrast to compact H2O maser sources 
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at 22 GHz, the 183 GHz maser emission is sp atially 
extended in high-mass (Orio n, Cer nicharo et al. Il994t 
W49 N, Gon zalez-Alfonso et al. 119951 SgrB2, Cernicharo 
et al. 120061) and even low-mass (HH7-11, Cernicharo 
et al. Il996l) star-forming regions. Due to maser amplifi- 
cation, the observed 183 GHz H20(3i3^22o) line bright- 
ness temperatures in these regions exceed those observed 
in 12 CO(J=1^0 ) by u p to a factor of 20 (e.g. for Orion, 
Cernicharo et al. 11994 Schulz et al. ll995|) . 

However, the H20(3i3^22o) line is strongly absorbed 
by the terrestrial atmosphere, rendering detection in 
nearby galaxies very difficult. The first successful extra- 
galactic detection of H?Q(3i3 — >2?n) has been reported 
in NGC 3079 (Humphreys et al. [2005) . Previous searches 
in low-z (U)LIR Gs ha ve only resulted in upper lim- 
its (Combes etal. I1997|) . Recently, the H 2 O(3i3^22 ) 
transition has also been detected towards Arp 220 (Cer- 
nicharo, Pardo & Weiss 2006, in prep., see Table 1). At 
high redshift, only the 1\\ — > 2o2 transition of para- 
H2O at 752.033 GHz (rest frame) has tentatively been 
detected in the z ~ 2.3 QS O IRA S F10214+4724 (Encre- 
naz et al. 119931 Casoli etal. ll994j) . 

In this paper, we report on a search for maser emis- 
sion from the 183 GHz water line in the strongly lensed, 
radio-loud z = 3.2 QSO MG 0751+2716 using the VLA 4 . 
Due to its strong magnification, MG 0751+2716 is the 
brightest CO source at high redshift that can be ob- 
served in the H20(3i3— >22o) transition with the VLA. 
Its L' co and Lfir (see Table 1) are comparable to lo- 
cal starburst galaxies like Arp 220. This suggests that 
this source is undergoing massive star-formation, which 
contributes significantly to its far-IR luminosity (e.g. 
Fig. 8 in Solomon & Vanden Bout 120051) . We use a 
standard concordance cosmology throughout, with Hq = 
71 k ms-^ Mpc- 1 , Q m = 0.27, and Q A = 0.73 (Spergel 
et al. l2003j) . 

2. OBSERVATIONS 

We observed the H 2 0(3i3^2 2 o) transition (f rest = 
183.3101GHz) towards MG 0751+2716 using the VLA 
in D configuration on 2004 June 19 and 27. At the tar- 
get z of 3.200, this transition is redshifted to 43.6453 GHz 
(6.87 mm). The total on-sky integration time amounts 
to 16 hr. Observations were performed in fast-switching 
mode using the nearby source 0745+241 for secondary 
amplitude and phase calibration. Observations were car- 
ried out under good weather conditions with 26 anten- 
nas. The phase stability in all runs was good (typically 
<25° rms for the longest baselines). For primary flux 
calibration, 3C286 was observed during each run. 

Two 50 MHz (corresponding to 344 km s" 1 at 
43.6 GHz) intermediate frequencies (IFs) were observed 
simultaneously in the so-called 'quasi-continuum' 
mode. One IF was centered on the line frequency at 
43.6351 GHz (closest possible tuning frequency, offset 
by -70kms~ 1 from the 12 CO(J=4^3) line center of 
Barvainis etal. I2002fl . and the second IF was centered 
at two frequencies symmetrically offset by 150 MHz 
(1030 kms" 1 , i.e. at 43.4851 GHz and 43.7851 GHz) 
from the line frequency to monitor the 7 mm continuum 

4 The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 



of MG 0751+2716 simultaneously. Despite the small off- 
set from the line center, the frequency/velocity coverage 
is well-matched to a line with a FWHM simila r to t he 
12 CO(J=4->3) line reported by Barvainis etal. l)2002j) . 

For data reduction and analysis, the AITS 5 package 
was used. The two continuum channels were concate- 
nated in the uv/visibility plane. The data were mapped 
using the CLEAN algorithm and 'natural' weighting; this 
results in a synthesized beam of 1.8"xl.7" (~13kpc at 
z = 3.2) at a major axis position angle of 78°. The final 
rms in both the line channel and the combined contin- 
uum channel is 60/iJy beam -1 . 

3. RESULTS 

In the top panel of Fig. 1, the final map of the 'ON' 
channel (line+continuum emission) is shown, while the 
middle panel shows the map of the two combined contin- 
uum ('OFF') channels, representing the continuum-only 
flux. The source is clearly detected in both maps. To de- 
rive the total flux density of the source, a 2-dimcnsional 
Gaussian was fitted to the detected source structure. 
Our best fit for MG 0751+2716 gives a Gaussian diam- 
eter of 2.3 "x 2.3" (source convolved with the synthe- 
sized beam); i.e. the emission appears only marginally 
resolved in our observations, as our beam is too large to 
recover the lens ima ge sub structure seen at higher res- 
olution (Carilli et al. I200"5|) . The integrated flux density 
in the 'ON' channel is 10.3 ± 0.2 mJy (peak: 6.12 ± 
0.06 mJy beam" 1 ), and the continuum flux in the com- 
bined 'OFF' channels is 10.3 ± 0.2 mJy (peak: 5.88 ± 
0.06 mJy beam" 1 ). From the VLA 42.2 GHz measure- 
ment of 13.2 mJy (Carilli et al. I2005D and using a spec- 
tral index of /3 = — 1.2 (Lehar et al. I1997J) . an extrapo- 
lated flux density of 12.7 mJy can be calculated, which 
is higher than derived from our observations. Therefore, 
continuum variability cannot be excluded. 

In the bottom panel of Fig. 1, the continuum- 
subtracted 'line' map is shown. This map was gener- 
ated by subtracting a CLEAN component model of the 
continuum emission (combined 'OFF' channel) from the 
visibility data of the 'ON' channel, which was then im- 
aged applying the CLEAN algorithm with the same pa- 
rameters used to create the continuum model. No clear 
evidence for H 2 0(3i3— >2 2 o) line emission is found. Thus, 
we set a 3cr upper limit of <Sh 2 o < 180 /iJy on the peak 
flux of H 2 0(3i3^2 2 o) emission. For MG 0751+2716, we 
thus derive an upper limit for the H20(3i3— >22o) line 
luminosity of L' H ^ < 0.6 x 10 9 K kms" 1 pc 2 (3tr limit, 
corrected f or gra vitational magnification, fi^ = 17, Bar- 
vainis et al. I2002L cf. Table [TJ . 

4. DISCUSSION 

4.1. Intrinsic line brightness temperatures 

Given that the 183 GHz water maser line most likely 
arises from extended r egion s associated with star forma- 
tion (Cernicharo etal. 1994), we can use our upper limit 
together w ith th e 12 CO(J=4— >3) line luminosity (Bar- 
vainis et al. 12002^1 to estimate how much of the molecular 
gas in MG 0751+2716 may be associated with massive 
star forming regions. Radiative transfer models based 
on the H20(3i3— >22o) emission line show that it is in- 
verted in warm and dense environments (Tki n > 40 K, 

5 www.aoc.nrao.edu/aips/ 
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nn 2 > 5 x 10 4 cm~ 3 ) which are typically found in star- 
forming regions, and that the resulting line brightness 
temperature critically depends on the underlying phys- 
ical conditions. It is known to range from thermal- 
ized emission up to line temperat ures o f ~10000K (Cer- 
nicharo et al. l!994t Combes et al. 119971 . As an example, 
observations of the 183 GHz water maser towards Orion 
IRc2 show that emission arises from narrow (few km s _1 ) 
line features with up to ~ 2000 K as well as a broad 
(~ 200 kms -1 ) emission plateau with ~ 500 K p e ak lin e 
brightness temperature (Cernicharo etal. 11994 11999ft . 
As an order of magnitude estimate, we assume in the 
following that the H 2 O(3i 3 ^2 20 ) line in MG 0751+2716 
arises from a similar environment as Orion with an intrin- 
sic brightness temperature of T b (H 2 0) ~ 500 K. Only the 
12 CO(J=4^3) line is observed so far in MG 0751+2716, 
therefore the brightness temperature of the CO emission 
is not well-constrained. However, the dust continuum 
temperature ofT^st ~ 40K in MG 0751+2716 (Barvai- 
nis & Ivison [2Q.Q2) is consistent with those observ ed in 
most other high-z QSOs (e.g. Beelen et al. 12006ft . In 
addition, the CO brightness temperature in well-studied 
high-z objects such as IRASF10214+4724 agrees within 
a fac tor of ~ 2 with the dust temperature (e.g. Downes 
etai ll995ft . By analogy, we here assume T b (CO) ~ 40K. 

4.2. Area filling factor 

The combined ratios of the brightness temperatures 
and line luminosities of H 2 and CO (T b (H 2 0) ~ 500 K, 
T b (CO) ~ 40K, L' H2Q < 0.6 x 10 9 K kms" 1 pc 2 , L' co = 
10.0 x 10 9 K kins -1 pc 2 ) can be used to estimate the rel- 
ative area filling factor (FF) of both molecules. By this 
means, our observed upper limit on the 183 GHz H 2 
line luminosity together with the 12 CO( J=4^3) line lu- 
minosity translates into an upper limit on the H 2 0/CO 
area filling factor of < 1%. 

To better understand the relevance of this filling factor 
limit derived from our observations, we now give an inde- 
pendent estimate for the filling factor based on the FIR 
luminosity and geometrical arguments. For this purpose, 
we assume in the following that all 183 GHz H 2 emis- 
sion in MG 0751+2716 is created in hot, star-forming 
cores like the one found in the central region of Orion. 
Assuming Lfir = 1.2 x 10 5 for t he central ID' of 
the Orion Nebula (Werner et al. 119761 corresponding to 
an area of 0.013 pc 2 at a distance of 450 pc), we find 
that 10 7 of these cores are needed to account for the far- 
IR luminosity in MG 0751+2716 (Table 1). Assuming a 
size of 0.013 pc 2 each, 10 7 of such cores would fill a disk 
with an equivalent radius of 200 pc. Taking the observed 
12 CO( J=4-^3) line FWHM and luminosity into account, 
and assuming T b (CO) ~ 40 K, for MG 0751+2716, we 
can calculate that the dust and CO are distributed over 
a region with an equivalent radius of 450 pc. Comparing 
the derived disk sizes, we thus find an expected H 2 0/CO 
area filling factor of ^20%. This estimate is much higher 
than what is actually observed. 

4.3. Comparison with other (U)LIRGs and NGC3079 

Our results for MG 0751+2716 are in line with observa- 
tions of the radio-quiet QSOs Mknl014 (PG 0157+001, 
z = 0.16) and VIIZw244 (PG0838+770, z = 0.13, see 
Table 1), the only two other infrared-luminous galaxies 



for which upper limits on th e 183 GHz H 2 line have 
been reported (Combes et al. 11997ft to date. Using the 
same assumptions on the intrinsic line brightness tem- 
peratures as for MG 0751+2716, the non-detections in 
these galaxies imply filling factors of the water line rela- 
tive to 12 CO(J=1^0) of < 5% and < 2% for Mknl014 
and VII Zw 244 respectively. As CO emission appears to 
be close to thermalized up to the 12 CO ( J=4^3) transi- 
tion in high-z QSOs (Weiss et al. l2005l) . these limits are 
directly comparable to MG 0751+2716. 

In Arp220 (z = 0.018), the H 2 O(3i 3 ^2 20 ) line was 
detected with a line luminosity of = 2.9 x 
10 8 Kkms _1 pc 2 (Cernicharo, Pardo & Weiss 2006, in 
prep.). With the CO luminosity given in Solomon et al. 
(1992a|, this translates to a H 2 0/CO filling factor of 
0.4%. Arp220 has a far-IR luminosity similar to that of 
MG 0751+2716, which might suggest a similar H 2 0/CO 
filling factor. This suggests that the H 2 0(3i3— >2 2 o) hne 
might be detectable in MG 0751+2716 if the sensitivity 
was increased only by a factor of few. 

For the recent first extragalactic H 2 0(3 13 ^2 20 ) detec- 
tion in the LINER NGC3079 (Humphreys etal. 2005), 
the H 2 0/CO luminosity ratio (using same brightness 
temperatures as before) translates into a H 2 0/CO area 
filling factor of only 0.05%. However, it is difficult 
to assess the relevance of the NGC 3079 results for 
MG 0751+2716, as the former is not a ULIRG. For 
lower-luminosity galaxies, it has been found that the 
HCN/CO luminosity ratio is signific antly lo wer (2.5-4% 
rather than 10-25%, Gao & Solomon 2004b); by analogy, 
also lower H 2 0/CO luminosity ratios may be expected 
in galaxies with FIR luminosities that are significantly 
lower than in a ULIRG. 

4.4. Possible interpretations of the non-detection 

Our observations suggest that the H 2 0/CO area fill- 
ing factor in MG 0751+2716 is significantly lower than 
expected from simple estimates based on the distribu- 
tion and temperature of the dust and CO. 

However, other reasons for the weakness of the wa- 
ter emission sho uld al so be considered: as discussed 
in Combes etal. (1997), the physical conditions to ef- 
ficiently pump the 183 GHz line to the assumed 500 K 
may not be present. There are strong differences in 
the conditions for H 2 line emission among Galactic 
star-forming regions. E.g., the H 2 0(3i3^2 2 o) maser 
emission in SgrB2(N/M) is much weaker than that in 
Orion IRc2. The extent of the maser emitting regions in 
both sources is of order 1 pc, but the observed bright- 
ness temperature in SgrB2 is by 1-2 orders of magni- 
tude lower , reac hing only few tens of Kelvin (Cernicharo 
etal. 11994 120061) . This is much lower than the 500 K 
we assumed for a typical H 2 emitting region. The 
183 GHz H 2 lines are thus only by a factor of 3 am- 
plified relative to CO in this source. For high densities 
and/or column densities, which are known to be present 
in the compact molecular distributions of luminous IR 
galaxies, the maser emission of the 183 GHz line may 
even be quenched. This would lead to thermalized emis- 
sion without maser amplification (Cernicharo et al. 11994 
Combes et al. Il99"7j) . For thermalized emission in an op- 
tically thick environment, the H 2 0(3i3^2 2 o) line is not 
expected to be stronger than those of other high density 
gas tracers such as HCN. In Arp 220, the HCN luminos- 
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ity is ^hcn = 9.4 x 10 8 K kms x pc 2 , i.e. by a factor 
of 3 brighter than H20(3i3— >22o) (Cernich aro, P ardo & 
Weiss 2006, in prep.; Gracia-Carpio et al. 12006ft . The 
H2O and HCN luminosities in Arp220 would thus even 
be consistent with Tb(H2 0) = 71, (HCN) and an area fill- 
ing factor that is smaller for H2O than for HCN. This 
suggests that only a small fraction of the de nse g as is 
giving rise to H2O emission. Carilli etal. (200]j) set 
a 3cr limit of L' HCN < 1-0 x 10 9 K kms~ 1 pc 2 on the 
HCN(J=2^1) line luminosity in MG 0751+2716, which 
is twice as high as our limit for L'^q. 

It also remains a possibility that the reason for 
our non-detection is due to limitations in our observ- 
ing mode. Due to the small fraction of gas giving 
rise to the H20(3i3— >22o) emission compared to CO, 
its linewidth could be much smaller than that of the 
12 CO(J=4^3) transition. This has r ecently been found 
in NGC3079 (Humphreys etal. 12005ft . where the differ- 
ence in linewidth is more than a factor of 10. As our 
observations had to be done using 50 MHz (344 kms" 1 ) 
channels, the emission from a narrow line would be di- 
luted over the full velocity range covered by that channel, 
rendering detection unlikely even at the achieved sensi- 
tivity. However, in Arp220, which is likely more similar 
to MG 0751+2716, the H 2 0(3i3^2 2 o) line has a width of 
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minous alternative to HCN to detect high-density gas in 
star-forming environments at high redshift. 
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TABLE 1 

EXTRAGALACTIC PARA-H20(3l3— >22o) EMISSION: FLUXES, LUMINOSITIES, AND AREA FILLING FACTORS. 
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D h 


^obs 


Sh 2 o 


J 1 a 


^H 2 


^FIR 


FF 


Ref. 






[Mpc] 


[GHz] 


[mjy] 


[lO'Kkms^pc 2 ] 


[10 9 K kms" 1 pc 2 ] 


[10 12 L S ] 






MG0751+2716 b 


3.200 


27940 


43.6453 


<0.18 


10.0 


< 0.6 


1.2 


<1% 


1,2,3 


Mkn 1014 


0.1631 


774 


157.6048 


<26 


7.8 


< 2.9 


2.2 


<5% 


4,5 


VII Zw 244 


0.1324 


616 


161.8775 


<26 


3.9 


< 0.7 


0.14 


<2% 


4,5 


NGO3079 


0.003723 


16 


182.6302 


550 


1.1 


0.0041 


0.021 


0.05% 


6,7 


Arp 220 


0.018126 


78 


180.0680 


170 


5.9 


0.29 


1.3 


0.4% 


8,9,10 



Ref erences. — [1] This w ork, [ 2] Barvainis et al. 1200211 . [3] Carilli et al. 120031 . [4] Combes et al. 119971). [51 Alloin etal. 119921) [6] K oda etal. 
1200211 . [7] Humphreys etal. 120051) . [8] Cernicharo, Pardo & Weiss (2006), in prep., [9] Solomon etal. I1992al) . [10] Solomon etal. 119971) . 



Note. — For MG0751+2716, the 12 CO(J=4^3) FWHM linewidth (Barvainis et al. 120021 corrected for the VLA b andpass) of 350 kms" 1 is 
utilized to derive L^q, whereas the 12 CO(J=1-*0) FWHM linewidths of 210 kms -1 and 80 kms -1 (Combes etal. 11997^1 arc assumed for Mkn 1014 
and VII Zw 244. 

a Derived as described by Solomon etal. I1992bl) : L' X [K km s~ *pc 2 ] = 3.25 X 10 7 X ! X X L>£ x (1 + z)~ 3 , where X is the molecule, / is the 

velocity— integrated line flux in Jy kms -1 , £>l is the luminosity distance in Mpc, and z^obs is the observed frequency in GHz. 
^This QSO is lensed by a factor of — 17 (Barvainis et al. l2002D . All given luminosities are corrected for lensing. 
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Fig. 1. — VLA observations of MG0751+2716 at a resolution of 1.8"xl.7" at a major axis position angle of 78°. The circle in the 
bottom left corner of each map represents the FWHM of the restoring CLEAN beam. Top: Map of the central 50 MHz (344 km s -1 ) 
'ON' channel at 43.6 GHz (183.3 GHz rest-frame), the frequency for which the H2O line is expected. Middle: Combined map of the two 
50 MHz 'OFF' (continuum only) channels symmetrically offset by 150 MHz (1030 km s — *) each from the H2O line frequency. The source is 
marginally resolved in both maps. Bottom: Continuum-subtracted ('ON— OFF') map. No clear evidence for H20(3i3^22o) line emission 
is found within the uncertainties of the observations. All three maps are shown with (-4, -3, 3, 4, 8, 16, 32, 48, 64, 80, 96)x60/iJy beam -1 
contours. 



